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Abstract 

A high-quality NIVIR solution structure is presented for protein hlVlcl-l(171-327) which comprises residues 171-327 of the 
human anti-apoptotic protein IVlcl-l (hMci-l). Since this construct contains the three Bcl-2 homology (BH) sequence motifs 
which participate in forming a binding site for inhibitors of hIVlcl-l, it is deemed to be crucial for structure-based design of 
novel anti-cancer drugs blocking the Mcll related anti-apoptotic pathway. While the coordinates of an NMR solution 
structure for a corresponding construct of the mouse homologue (mMcl-1) are publicly available, our structure is the first 
atomic resolution structure reported for the 'apo form' of the human protein. Comparison of the two structures reveals that 
hIVlcl-l (171 -327) exhibits a somewhat wider ligand/inhibitor binding groove as well as a different charge distribution within 
the BH3 binding groove. These findings strongly suggest that the availability of the human structure is of critical importance 
to support future design of cancer drugs. 
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Introduction 

The malfunctioning of cellular apoptosis [1] is a major hallmark 
of cancer. The regulation of apoptosis depends on the family of 
Bcl-2 proteins which contain one or several Bcl-2 homology (BH) 
sequence motifs. Based on their function and the similarity of their 
respective BH sequence motifs, these proteins can be grouped into 
three classes [2], [3]: (i) multi-domain pro-apoptotic proteins such 
as Bax and Bak, (ii) anti-apoptotic (i.e., pro-survival) proteins such 
as Mcl-1, Bcl-1, Bc1-xl, Bcl-w and Bfl-l/Al, all of which exhibit a 
similar architecture as Bax and Bak, and (iii) several pro-apoptotic 
proteins comprising only a single BH3 sequence motif such as Bid, 
Bad, Bim, Puma, Noxa, Hrk, Bmf, and Nbk/Bik ('BH3-only' 
proteins). The BH3 motif of class (iii) proteins forms an 
amphipathic ot-helix which interacts specifically with a hydropho- 
bic pocket formed in both pro-apoptotic class (i), and anti- 
apoptotic class (ii) proteins with participation of their respective 
BH motifs [2], [3]. Inhibition of the resulting protein-protein 
complex formation offers a promising strategy to treat cancer. For 
example, the small molecule Bcl-2 antagonist ABT-737 [4] inhibits 
anti-apoptotic class (ii) proteins Bc1-xl, Bcl-w and Bcl-1, and a 
congener [5] that can be orally administered is currently in clinical 
trials. 

The anti-apoptotic, pro-survival 350-residue protein Mcl-1 
('myeloid cell leukemia- 1') [2] is primarily anchored in the outer 
mitochondrial membrane by a C-terminal trans-membrane 
domain and contains three BH sequence domains: BH3 (residues 



209-223), BHl (252-272) and BH2 (residues 304-319) [2]. Mcl-1 
inhibits death receptor-induced apoptosis by selectively binding to 
truncated Bid (tBid) [6] and can sequester endogenous Bak to 
block Bak-mediated cell death. Moreover, Mcl-1 interacts with 
several BH3-only proteins (Bim, Bid and Puma, Noxa and Bak). 
Hence, Mcl-1 plays an early role in response to signals directing 
either cell survival or cell death [2] and has been shown to be up- 
regulated in numerous malignant tumors. Approaches abrogating 
the Mcl-1 's anti-aptototic function either by reducing its abun- 
dance or by inactivating its functional BH3-binding groove show 
great promise for the cancer treatment [2] , [4] , [6] , [7] . Here we 
present the high-quality NMR solution structure of polypeptide 
segment 171-327 of human Mcl-1 (hMcl-1) which comprises the 
three BH motifs deemed to be crucial for structure based drug 
design. 

Results and Discussion 

A high-quality NMR structure of hMcl-l(171-327) was 
obtained (Table 1) and the coordinates were deposited in the 
PDB [8] (accession code 2mhs). The structure comprises seven ot- 
helices al-ot7 (residues 173-191, 204-235, 240-253, 262-280, 
284-301, 303-308 and 31 1-319) arranged to form the character- 
istic 'Bcl-2 core' structure [9] (Figure 1). The helices are locally 
and globally well-defined, while the C-terminus (residues 320-327) 
and the loops connecting, respectively, helices al and a2, helices 
a3 and a4, and helices a4 and a5 are flexibly disordered. The 
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Figure 1. NMR structure of hMcl-1(171-327). (A) Backbone of the 
20 CYANA conformers representing the solution structure of hMcl- 
1(171-327) after superposition of bacl<bone N, and C atoms of the cx- 
helices for minimal rmsd. The three BH sequence motifs are colored in 
green (BH3), red (BH1) and blue (BH2), respectively. (B) Ribbon drawing 
of the lowest energy conformer of hMcl-l (171-327). a-helices a1-0(7 are 
labeled and colored differently, and the N- and C-termini are labeled as 
"W " and "C ". The figures were generated using the programs 
MOLiVlOL [36] and PYMOL [37]. 
doi:10.1371/journal.pone.0096521.g001 

central helix a4 is surrounded by the other six helices, with otl, a2, 
a3 and 0(5 packed around one side, and 0(6 and a? packed against 
its N-terminus. Helices 0(2, Q(3, 0(4 and 0t7 participate in forming 
the BH3 binding groove. The electrostatic protein surface 
potential is positive at both ends of the BH3 binding groove 
(due to the presence of Arg 233, Lys 234, Arg 248 and Arg 263) 
and negative at the side of helix a3 side (due to Asp 256) (Figure 2). 
This shows that the charge distribution in the BH3 binding groove 
of h]V[cl-l(171-327) differs distinctly from other anti-apoptotic 
proteins [10]. 

Including our hMcl-l(171-327) structure, twenty atomic 
resolution structures containing different Mcl-l constructs are 
currently deposited in the PDB. In addition to the two 'apo' 
proteins h]VIcl-l(171-327) and mouse m]VIcl-l(152-308) [10] 
[PDB accession code Iwsx, 89% sequence identity with the 
human protein], the structures for nineteen protein-ligand 
complexes were deposited (Table 2) [9], [11-18]. Clearly, the 
large number of available structures reflects the outstanding 
interest in Mcl- 1 as a target for the development of new cancer 
drugs. Superposition of the ot-hehces reveals, as expected, close 
structural similarity for all Mcl-l proteins structures (Figure 3): the 



root mean square deviation (rmsd) values range from 1.05 to 
1.54 A relative to hMcll-l(171-327) (Table 2). However, com- 
parison of the two apo protein structures of hMcl-l(171-327) and 
m]VIcl-l(l 52-308) with the complex structures shows that the 
binding pocket is widened upon complex formation (Table 2): the 
distances between the C'^-atoms of residues His 224 in helix ot2 
(His 205 in mMcl-1) and His 252 (His 233 in mMcl-1) at the C- 
terminus of helix 0(3 are, respectively, ~ 1 6 A and ~ 1 4 A in hlVIcl- 
1(171-327) and mMcl-l(152-308), and -18-21 A in the 
complexes. 

The fact that the human apo protein exhibits a somewhat wider 
binding groove than the mouse homologue (Table 2) can be, at 
least partially, ascribed to the side chain of Leu 246 in the human 
protein which is not buried as deeply as the corresponding Phe 
side chain in the mouse protein. Furthermore, when comparing 
the human and the mouse protein, differences are observed for the 
charge distributions in the BH3-binding groove (Figure 2): the 
human protein is negatively charged on the side of helix 0(3, while 
the corresponding surface of mouse protein is positively charged. 
This difference arises from Ser 255 corresponding to Lys 236 in 
the mouse protein. Remarkably, hMcl-l(171-327) is structurally 
more similar to the h]VIcll(l 7 l-327)-hBim BH3 complex (Figure 3) 
than to apo mMcl-1 (152-308) (Table 2). 

Taken together, structural comparisons show that, in spite of the 
89% sequence identity between human and mouse protein, the 
availability of the human hMcl-l(171-327) structure can be 
expected to be of critical importance for supporting future design 
of cancer drugs. 

Materials and Methods 

NMR Sample Preparation 

Preliminary studies showed that hMcl-l(171-327) (UniProtKB/ 
Swiss-Prot ID QP7820/MCL1_HUMAN) is not stable in solution. 
However, the mutant Cys 286 Ser is stable for several weeks at 
concentrations ~0.7 mM, and both wild- type and mutant bind 
the Bim-BH3 peptide with the same affinity (K^ ~ 60 pM) in a 
Biacore assay. Hence, we solved the NMR structure of hMcl- 
1(171-327) Cys 286 Ser referred to as hMcl-l(17 1-327) in this 
publication. 

hMcl-l(l 7 1-327) was cloned, expressed, refolded and purified 
following standard protocols to produce a uniformly '^C, '^N- 
labeled protein sample [19]. Briefly, the gene was cloned into a 
pET21d (Novagen) derivative yielding plasmid pSR482-21.1. The 
resulting construct contains seven nonnative residues at the C- 




A B 

Figure 2. Electrostatic surface potentials. (A) For human hMcl-1(171-327) in the orientation shown in Figure 1 (left) and after rotation by 180° 
about the vertical axis (right). Surface colors (blue for positively charged; red for negatively charged) indicated the electrostatic potential calculated 
by using PYMOL [37] and its default vacuum electrostatics protocol. (B) Same as in (A) but for mouse miVlcl-1 (152-308). 
doi:1 0.1 371/journal.pone.0096521 .g002 
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Table 1. Statistics of hMcl-l{171-327) NMR Structure. 





Completeness of stereo-specific assignments [%]^ 


"CH2 of Gly 


55 (6/11) 


l'CH2 


38 (27/71) 


Val and Leu methyl groups 


100 (27/27) 


Conformationally restricting distance constraints 


Intra-residue [i = j] 


1052 


Sequential [|i-j| =1] 


1062 


Medium range [1 < |i-j| < 5] 


1197 


Long range [[i-j| ^ 5] 


1058 


Total 


4369 


Didedral angle constraints 


cp 113 




113 


Number of constraints per residue (170-327) 


29.1 


Number of long range constraints per residue (170-327) 


6.7 


CYANA target function [A^] 


0.88±0.12 


Number of distance constraints violations per CYANA conformer 


0.2-0.5 A 


0 


> 0.5 A 


0 


Number of dihedral-angle constraint violations per CYANA 
conformer 


> 5 


0 


Average rmsd to the mean CNS coordinates [A] 


A-helices,'^ backbone heavy atoms N, C^ C 


0.42 ±0.05 


A-helices,'^ all heavy atoms 


0.88±0.07 


Residues 172-312, backbone heavy atoms N, C", C 


0.65±0.13 


all residues, all heavy atoms 


1.05±0.10 


PROCHECK [38] G-factors raw score 


((p and \|//all dihedral angles)^ 


0.34/0.22 


PROCHECK [38] G-factor Z - score 


((p and \|//all dihedral angles)*^ 


1.65/1.30 


MOLPROBITY[39] clash score (raw/Z - score)' 


20.88/- 2.06 


AutoQF R/P/F/DP scores [40] (%) 


96/97/96/81 


Ramachandran plot summary'^ 


most favorable regions 


92.7 


additionally allowed regions 


7.3 


generously allowed regions 


0.0 


disallowed regions 


0.0 



*^Related to pairs with non-degenerate chemical shift. 

''Regular secondary element: a-he!ical residues 173-191, 204-235, 240-253, 

262-280, 284-301, 303-308 and 311-319. 

'Ordered residues: 172-192,194-198, 204-235, 238-255, 262-321 with dihedral 
angle order parameters S((p) and S(\|;) > 0.90. Z-scores were computed relative 
to corresponding structure quality measures for high resolution X-ray crystal 
structures [42]. 

doi:l 0.1 371/journal.pone.0096521 .tool 

terminus (LHHHHHH) to facilitate protein purification. Esche- 
richia Coll BL2 1 (DE3) pMGK cells, a codon enhanced strain, were 
transformed with pMcl 1 -2 1 . 1 , and cultured in MJ9 minimal 
medium [20] containing ('^NH4)2S04 and f/-'^C-gl ucose as sole 
nitrogen and carbon sources. Double-washed inclusion bodies 
containing hMcl-l(171-327) were solubilized in 8 M buffered 
guanidine hydrochloride (VWR) containing 5 mM DTT, slowly 



C 




A B C 



Figure 3. Superposition of selected Mcl-1 structures. (A) 

Structures of hMcl-1(171-327) (green) and mMcl-1(1 52-308) (cyan, 
PDB accession code Iwsx) after superposition of the backbone N, C" 
and C atoms of the ot-helices for minimal rmsd. (B) Ribbon drawing 
(zoomed into (A)) showing the different binding groove widths of 
human (green) and mouse (cyan) protein. The distances between the 
Ca-atoms of residues His 224 in helix a.2 (His 205 in mMcl-l ) and His 252 
(His 233 in mMcl-1 ) at the C-terminus of helix cx3 are highlighted: ~ 1 6 A 
in hMcl-1 (1 71-327) and ~ 14 A mMcl-1 (1 52-308) (C) Superposition as in 
(A) of hMcl-1(171-327) (green) and mMcl-1 (152-308) (cyan, Iwsx), and 
six selected Mcl-I complex structures (see also Table 2): human Mcl-1 
complexed with Bim BH3 (magenta, 2nla); chimeric rat-human rMcl- 
1(171-208)hMcl-1 (209-327) complexed with mouse mNoxaB BH3 
(yellow, 2rod); mouse mMcl-1 (152-308) complexed with mouse NoxaA 
BH3 (pink, 2roc); mouse mMcl-1 (1 52-308) complexed with mouse Puma 
BH3 (grey, 2jm6); mouse mMcl-1 (152-308) complexed with mouse 
NoxaB BH3 (purple, 2rod); chimeric rat-human mMcl-1 (171 -208)hMcl- 
1(209-327) complexed with human Bim BH3 (orange, 2nl9); chimeric 
rat-human mMcl-1 (171-208)hMcl-1 (209-327) complexed with human 
Bim BH3(L62A, F68A) (light green, 3d7v).The figures were prepared with 
the programs MOLMOL [36] and PYMOL [37]. 
doi:1 0.1 371/journal.pone.0096521 .g003 

diluted into nine volumes of 20 mM Tris-HCl, 250 mM NaCl, 
0.5 M urea, 10% glycerol, pH 7.4, and refolded within 72 hours 
at 4°C. hMcl-l(171-327) was purified using a Talon affinity resin 
(Clontech) applied to a HiTrap SP High Performance column (GE 
Healthcare). The final yield of purified U-^'^G, ^^N protein (> 98% 
homogeneous by SDS-PAGE; 20.3 kDa by MALDI-TOF mass 
spectrometry) was ~25 mg/L. In addition, a fZ-'^N and 5% 
biosynthetically directed fractionally ''^C-labeled sample [21] was 
generated for stereo-specific assignment of isopropyl methyl 
groups. NMR samples were prepared at 0.7 mM protein 
concentration. An isotropic overall rotational correlation time of 
~10 ns was inferred from ' 'N spin relaxation times indicating that 
hMcl-l(171-327) is monomeric in solution. 

NMR spectroscopy 

NMR spectra were recorded at 25°C. Five G-matrix Fourier 
transform (GFT) NMR experiments [2 2], [2 3] and a simultaneous 
3D ''N/'''C''"P'''"'V"C""""""'^-resolved NOESY [24], [25] spec- 
trum (mixing time 60 ms; measurement time: 48 hours) were 
acquired on a Varian INOVA 750 MHz spectrometer equipped 
with a conventional probe. 2D constant-time [''^C,'Fr|-HSQC 
spectra (18 hours) were recorded for the 5% biosynthetically 
directed fractionally ''^C-labeled sample on a Varian INOVA 
600 MHz spectrometer equipped with a cryogenic probe as was 
described [21], [26]. Spectra were processed and analyzed using 
the programs NMRPipe [27] and XEASY [28]. 

Sequence specific backbone (H^, H", N, C") and H'^/C'' 
resonance assignments were obtained by using (4,3)D HNNC^'^C 
(63 hours)/(4,3)D C°'''C°'(CO)NHN (62 hours), and (4,3)D H^- 
C'"'' (CO)NHN (69 hours) [23] along with the program AUTO- 
ASSIGN [29]. More peripheral side chain chemical shifts were 
assigned with aliphatic (4,3)D HCCH (87 hours) [23] and 3D 

l,5J^/l;^^aliphatkyl:^^ar„matk,_^g^^^^gJ ['H, 'H] -NOESY [24], [25]. 
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Table 2. Rmsd values for comparison of the NMR structure of hMcl-l(171-327) with the structures of mouse mMcl-l (152-308) and 
Mcl-l complexes.^ 



172-193, 203- 



Mcl-1 structures 


321" 


a1 -a?" 


209-321"" 


a2-a4'' 


02-07'' 


dCA224,252'' 


mMcl-l" 


1.60 ±0.09 


1.52±0.06 


1.61±0.10 


1.21 ±0.07 


1.53 ±0.06 


13.2-14.6 


hMcl-l-hBim' 


1.76±0.10 


1.41 ±0.06 


1.87±0.11 


1.52±0.09 


1.53 ±0.09 


19.9 


rMd-l-hMcl-l-mNoxaBS 


1.46±0.12 


1.05 ±0.08 


1.53±0.14 


1.00±0.10 


1.08±0.09 


19.9 


mMcl-l-mNoxaA ^ 


1.52±0.11 


1.16±0.07 


1.59±0.13 


i.ii±c.n 


1.18±0.08 


18.8-20.2 


mMcl-1-mPuma' 


1.57±0.09 


1.30±0.08 


1.59±0.11 


1.24±0.11 


1.30±0.09 


18.3-19.9 


mMcl-1-rrNoxaB' 


1.46 ±0.09 


1.13±0.05 


1.53±0.11 


1.12±0.08 


1.18±0.06 


18.3-19.6 


rMd-1-hMcl-l-hBim'' 


1.75±0.09 


1.44 ±0.06 


1.86±0.11 


1.58±0.09 


1.56±0.07 


19.9 


rMcl-1-hMcl-l-hBim(L62A, F78A) ' 


1.80 ±0.09 


1.44±0.06 


1.90±0.11 


1.53±0.08 


1.55±0.07 


19.8 


hMcM-hBidBH3-' 


1.84±0.15 


1.54 ±0.09 


1.57±0.20 


1.30±0.15 


1.30±0.12 


20.8-21.4 


hMcM-hBIMBHB" 


1.38±0.16 


1.06±0.07 


1.43±0.18 


0.93±0.10 


1.09±0.09 


19.4 


hMcl1-BimL12Y' 


1.75±0.13 


1.50±0.08 


1.85±0.15 


1.48±0.12 


1.59±0.09 


20.2 


hMcll-BimBH3 2dA"' 


1.75±0.13 


1.50±0.08 


1.84±0.14 


1.46±0.12 


1.57±0.09 


19.7 


hMcl1-BimBH3 F4aE" 


1.75±0.13 


1.47 ±0.08 


1.84±0.15 


1.44±0.08 


1.55±0.08 


19.9 


hMcll -B7° 


1.73±0.14 


1.46±0.08 


1.83±0.11 


1.42±0.11 


1.54±0.08 


19.6-19.9 


hMcM-hMcMBHB" 


1.47±0.16 


1.21 ±0.09 


1.48±0.18 


1.08±0.15 


1.24±0.11 


20.3 


hMcll-BaxBHS" 


1.69±0.12 


1.46±0.08 


1.79±0.13 


1.43±0.12 


1.54±0.09 


19.4 


mMcll-NoxaBHa' 


1.43±0.16 


1.20±0.10 


1.48±0.17 


1.10±0.13 


1.23±0.11 


19.1 


hMcll -compound53^ 


1.45±0.12 


1.22±0.08 


1.48±0.14 


1.19±0.10 


1.28±0.08 


18.7 


hMcll -compound60^ 


1.45±0.13 


1.22 ±0.08 


1.47±0.12 


1.12±0.09 


1.25 ±0.08 


17.9-19.6 


hMcM-BH3" 


1.51±0.16 


1.22 ±0.09 


1.57±0.18 


1.06±0.12 


1.27±0.10 


20.3 



Average pairwise rmsd values {A) were calculated for backbone heavy atoms N, C", and C between the 20 conformers of Mcl-1 (171 -327) and corresponding 
polypeptide segments in the other structures. The distances dCA (in A) between the C"-atoms of residues His 224 in helix a2 (His 205 in mMcl-1) and His 252 (His 233 in 
mMcl-1) at the C-terminus of helix 0(3 are provided as a measure for the width of the BH3 binding groove. 

"^Residue numbers are for hMcl-l(171-327); residues 194-202 were excluded since one structure (2nl9^) does not contain the corresponding residues; residues 172-193 
and 203-321 correspond to residues 153-174 and 184-302 In mMcl-1, and residues 209-321 correspond to residues 190-302 in mMcl-1. 

^Helices al-a7 In hMcl-1 comprise residues 173-191, 204-235, 240-253, 262-280, 284-301, 303-308 and 31 1-319; the corresponding residues in mMcl-1 are: 155-172, 
185-216, 221-234, 243-261, 265-282, 284-289 and 292-300. 

^Helices a2-a7 In hMcl-1 and residues 204-208 (numbers In hMcl-1) were excluded 

"Mouse mMcl-1 (152-308), PDB accession code Iwsx (the mean NMR coordinates were used) [10]. 

'Human hMcl-1 complexed with human hBim BH3, 2pqk [11]. 

^Chimiric rat-human rMcl-l(l 71 -208)hMcl-l (209-327) complexed with mouse mNoxaB BH3, 2nla [9]. 

'^Mouse mMcl-1 complexed with mouse mNoxaA BH3, 2rod [12]. 

'Mouse mMcl-1 complexed with mouse mPuma BH3, 2roc [12]. 

^Mouse mMcl-1 complexed with mouse mNoxaB BH3, 2jm6 [9]. 

^Chimlrlc rat-human Mcl-1 complexed with human hBim BH3, 2nl9 [9]. 

'Chlmlric rat-human Mcl-1 complexed with human hBim (L62A, F68A), 3d7v [13]. 

^Human hMcll complexed with human Bid BH3, 2kbw [15]. 

^Human hMcl-1 complexed with human Bim BH3 mutant l2dY, 3kj0 [11]. 

'Human hMcl-1 complexed with human BimL12Y, 3io9 [16]. 

""Human hMcll complexed with human Bim BH3 mutant l2dA, 3kj1 [11]. 

"Human hMcll complexed with human Bim BH3 mutant F4aE, 3kj2 [11]. 

°Human hMcl-1 complexed with Mcll specific selected peptide B7, 3kz0 [41]. 

''Human hMcl-1 complexed with human Mcll BH3, 3mk8 [17]. 

^Human hMcll complexed with human Bax BH3, 3pkl. 

^Mouse mMcl-1 complexed with mouse Noxa BH3, 4g35 [18]. 

^Human hMct-1 complexed with 6-chloro-3-[3-(4-chloro-3,5-dlmethylphenoxy)propyl]-1 H-indole-2-carboxylic acid, 4hw2 [14]. 
^Human hMcl-1 complexed with 6-chloro-3-[3-(4-chloro-3,5-dimethylphenoxy)propyl]-lH-indole-2-carboxylic acid, 4hw3 [14]. 
"Human hMcl-1 complexed with human Mcll BH3, 4hw4 [14]. 
doi:1 0.1 371 /journal.pone.0096521 .t002 



Overall, assignments were obtained for 96% of backbone and 
IjjP/ISqP resonances and for 93% of the side chain resonances 
which are assignable with the NMR experiments listed above 
(excluding the N-terminal NHa""", Pro ' 'N, ''^C preceding prolyl 
residues, Lys NHj""", Arg NH2, OH, side chain ' C and aromatic 
^^C^). Furthermore, 100%/100% of Val and Leu isopropyl 
moieties with non-degenerate proton chemical shifts were stereo- 
specificaUy assigned (Table 1). Chemical shifts were deposited in 



the BioMagResBank [30] (accession code 19654). 'H-'H upper 
distance limit constraints for structure calculations were obtained 
from NOESY (Table 1). In addition, backbone dihedral angle 
constraints were derived from chemical shifts using the program 
TALOS [31] for residues located in weU-defmed secondary 
structure elements (Table 1). The programs CYANA [32], [33] 
and AUTOSTRUCTURE [34] were used in paraUel to assign 
long-range NOEs [24]. The final structure calculations were 
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performed using CYANA followed by explicit water bath 
refinement using the program CNS [35]. 
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